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Abstract 

The solar system is moving through the partially 
ionized local interstellar cloud (LIC). The ionized 
matter of the LIC interacts with the expanding solar 
wind forming the heliospheric interface. The neutral 
component (interstellar atoms) penetrates through 
the heliospheric interface into the heliosphere, where 
it is measured directly "in situ" as pick-up ions and 
neutral atoms (and as anomalous cosmic rays) or in- 
directly through resonant scattering of solar Ly a. 
When crossing the heliospheric interface, interstellar 
atoms interact with the plasma component through 
charge exchange. This interaction leads to changes of 
both neutral gas and plasma properties. The helio- 
spheric interface is also the source of radio emissions 
which have been detected by the Voyager since 1983. 
In this paper, we have used a kinetic model of the 
flow of the interstellar atoms with updated values of 
velocity, temperature, and density of the circumsolar 
interstellar hydrogen and calculated how all quanti- 
ties which are directly associated to the observations 
vary as a function of the interstellar proton number 
density n P] Lic- These quantities are the degree of fil- 
tration, the temperature and the velocity of the inter- 
stellar H atoms in the inner heliosphere, the distances 
to the bow shock (BS), heliopause, and termination 
shock, and the plasma frequencies in the LIC, at the 
BS and in the maximum compression region around 
the heliosphere which constitutes the "barrier" for ra- 
dio waves formed in the interstellar medium. Com- 
paring the model results with recent pickup ion data, 
Ly a measurements, and low- frequencies radio emis- 
sions, we have searched for a number density of pro- 
tons in the local interstellar cloud compatible with all 
observations. 

We find it difficult in the frame of this model with- 
out interstellar magnetic field to reconcile the distance 
to the shock and heliopause deduced from the time de- 
lay of the radio emissions with other diagnostics and 



discuss possible explanations for these discrepancies, 
as the existence of an additional interstellar magnetic 
pressure (2.1 /j,G < B < 4 /j,G for a perpendicular 
magnetic field). We also conclude that on the basis 
of this model the most likely value for the proton den- 
sity in the local interstellar cloud is in the range 0.04 
cm~ 3 < n P! Lic < 0.07 cm -3 . 



1. Introduction 

Our solar system is moving through a partially ion- 
ized interstellar cloud. The ionized fraction of this lo- 
cal interstellar cloud (LIC) interacts with the expand- 
ing solar wind and forms the LIC-solar wind (SW) 
interface (or heliospheric interface). The characteri- 
zation of this interface is a timely major objective in 
astrophysics and space plasma physics. The interest 
to the construction of the LIC/SW interaction mod- 
els is increasing at the present time [Ripken and Fahr, 
1983; Baranov and Malama, 1993; Zank et al., 1996; 
Linde et al., 1998; Pogorelov and Matsuda, 1998]. The 
choice of an adequate model of the interface depends 
on the parameters of the LIC. Some of these param- 
eters, as the Sun/LIC relative velocity, or the LIC 
temperature are now well constrained [Witte et al., 
1993; Lallement and Bertin, 1992; Linsky et al., 1993; 
Lallement et al, 1995; see, also, Frisch, 1995], but un- 
fortunately there are no direct ways to measure the 
circumsolar interstellar electron (or proton) density, 
nor the local interstellar magnetic field, while these 
two parameters govern the structure and the size of 
our heliosphere. There have been measurements of 
the average electron density in the LIC toward nearby 
stars. However, resulting densities range from 0.05 (- 
0.04, +0.14) cm- 3 to up to 0.3 (-0.14, +0.3) cm" 3 
depending on the ions used for the diagnostics or on 
which linc-of-sight is probed [e.g., Lallement and Fer- 
let, 1997]. The most precise, temperature indepen- 
dent value is 0.11 cm -3 toward the star Capella [ Wood 
and Linsky, 1997]. In addition, what is measured is 
always averaged over large distances, while the ion- 
ization degree in the local interstellar medium is very 
likely highly variable and out of ionization equilibrium 
[e.g., Vallerga, 1998]. Therefore there is a need for 
indirect observations which can bring stringent con- 
straints on the plasma density and on the shape and 
size of the interface. Such constraints should help to 
predict when the two Voyager spacecraft will cross 
the interface and whether or not they will be able to 
perform and transmit direct observations. 
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Among the various types of heliospheric interface 
diagnostics, there are measurements of the pick-up 
ions. Pick-up ions are formed when interstellar neu- 
trals, having penetrated into the heliosphere, become 
ionized by charge exchange with the solar wind ions or 
photoionized. The newly created ions are then con- 
verted away from the Sun by the solar wind. The 
detection of He + and He ++ provides a new determi- 
nation of the neutral helium flow properties [Gloeck- 
ler et ai, 1997], which can be compared with the 
direct detection of the neutral helium [Witte et ai, 
1993, 1996]. It is interesting to note that the dif- 
ferent determinations of the helium density are now 
in rather good agreement. With the Solar Wind Ion 
Composition Spectrometer (SWICS) instrument on 
board Ulysses, Gloeckler et al. [1997] found n(Hel) 
= 0.0153 ± 0.0018 cm' 3 with an uncertainty of the 
order of only 8%, compatible with Active Magneto- 
spheric Particle Tracer Explorers (AMPTE) results 
of Moebius [1996] and Ulysses Interstellar Neutral- 
Gas instrument (GAS) results of Witte et al. [1996]. 
As for hydrogen, the first successful detection of H + 
pick-up ions has been done with the SWICS instru- 
ment on board the Ulysses spacecraft [Gloeckler et 
al., 1993]. From the H + fluxes, one can infer a new 
value for the neutral hydrogen flux in the inner helio- 
sphere. Indeed, Gloeckler et al. [1997] used the pick- 
up ions data to infer the interstellar hydrogen atom 
number density in the outer heliosphere and found 
n(HI) = 0.115 ± 0.025 cm' 3 with an uncertainty of 
20%. Note that in all these works the H and He num- 
ber densities are obtained using the classical so-called 
"hot" model of interstellar neutrals flow in the helio- 
sphere. The H number density obtained from pickup 
ion measurements is close to the lower range of the 
interval derived from optical resonance data revised 
by Quemerais et al. [1996], i.e., 0.11-0.17 cm' 3 . 

The neutral H density in the inner heliosphere is 
dependent on the perturbations the neutral H suf- 
fers at the heliospheric interface. Since interstellar 
helium is not depleted in the heliospheric interface 
region (n(Hc) heliospheric n(He) interstellar), it is 
possible to compare the properties of the neutral hy- 
drogen and helium flows. Taking into account the 
newly derived neutral hydrogen to neutral helium ra- 
tio in the LIC [Dupuis et al, 1995] 



n L i C (HI) 
riLic(Hel) 



= 14 ± 1 



proxy for the filtration ratio of H as a function of the 
plasma density. The proxy was derived from Bara- 
nov and Malama [1993] model filtration ratios. It 
was found that if n(H) in the inner heliosphere is of 
the order of 0.15 (0.10, respectively) cm' 3 , then the 
electron density in the circumsolar medium is close to 
0.05 (0.11, respectively) cm' 3 . Gloeckler et al. [1997], 
using the new SWICS pick-up ions results and an in- 
terstellar HI/Hel ratio of 13 ± 1 (the average value of 
the ratio toward the nearby white dwarfs), concluded 
that nLic(HI) = 0.2 ± 0.03 cm -3 , which corresponds 
to a filtration factor (the ratio of atom number den- 
sity in the outer heliosphere to atom number density 
in the LIC) 



n T s(HI) 
X = 



^-^- = 0.58 ±0.15 
nLic(HI) 



(1) 



Lallement [1996] related the neutral H density in the 
heliosphere to the plasma density in the LIC, using a 



where subscript TS is termination shock. Then, on 
the basis of estimates of the charge-exchange pro- 
cesses, they obtained for the interstellar proton (or 
electron) number density n Pi Lic = 0.04± 0.017 cm' 3 . 

For a given plasma density, there is a nonneg- 
ligiblc influence of the neutral density in the LIC 
on the filtration ratio. The filtration ratios used 
by Lallement [1996] were taken from a model with 
^lic(HI) = 0.14 cm' 3 [Baranov and Malama, 1993], 
introducing some unconsistency in the method. In 
what follows, we will make use of the more appropri- 
ate value nLic(HI) = 0.2 cm' 3 , which is based on the 
well-measured neutral helium density and the inter- 
stellar ratio measured with the Extreme Ultraviolet 
Explorer (EUVE). 

The goal of this paper is to determine a range for 
the interstellar proton number density n p ,Lic that 
is compatible with all observations, using the two- 
shock heliospheric interface model of Baranov and 
Malama [1993, 1995, 1996] for the updated value 
n LIC (HI) = 0.2 cm' 3 . The observations we will con- 
sider are, in addition to the pick-up ions quoted above, 
the temperature and velocity of the neutral H flow in 
the inner heliosphere and the Voyager radio emissions. 
Recently, Linsky and Wood [1996] have detected the 
heated and decelerated gas from the so-called H wall 
corresponding to the compressed region between the 
bow shock (BS) and the heliopausc (HP), in absorp- 
tion toward the star a Centauri. Gayley et al. [1997] 
have compared the observed absorption with the the- 
oretical absorption for three different models. One of 
these models is a two-shock model, and the two other 
correspond to the "subsonic" case; That is, they have 
modified the equation of state of the gas to simulate 
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the effect of an interstellar magnetic field. These au- 
thors conclude that the H wall absorption favors the 
"subsonic case." We have not included such diagnostic 
in our study for the following reasons: in our "super- 
sonic" case, the simulations show that it is hard to dis- 
tinguish H walls built up for n p lic = 0.04 cm~ 3 or for 
0.2 cm~ 3 . As a matter of fact, if n Pi Lic increases, the 
gas is more heated and compressed, but the thickness 
of the H wall is reduced. Also, the precision required 
to model the differences between the theoretical ab- 
sorptions, namely small differences of the order of a 
few kilometers per second at the bottom of the lines, 
is of the order of the differences between the model 
results from different groups for the same parameters 
in the supersonic case [see Williams et ah, 1997, Ap- 
pendix B]. A larger difference may exist between the 
"supersonic" and the "subsonic" cases, large enough 
to favor the subsonic case, as argued by Gayley et 
al. Our approach is to use other independent di- 
agnostics, namely all heliospheric data, in order to 
constrain the requirements for additional physics. 




Figure 1. Qualitative picture of the solar wind interaction with 
the local interstellar medium (LISM). Here BS is the bow shock, 
HP is the heliopausc, and TS is the termination shock. 



The model is a self-consistent gasdynamic-kinetic 
model of the solar wind interaction with the local in- 
terstellar medium, which takes into account the mu- 
tual influence of the plasma component of the LIC 



and interstellar H atoms in the approximation of ax- 
ial symmetry (Figure 1). 

Then we compare the range of « p ,lic values de- 
rived from the pick-up measurements with the range 
that is compatible with observations of the backscat- 
tered solar Ly a emission as well as from the interpre- 
tation of the 2-3 kHz emission recorded by the Voy- 
ager radio instruments. As a tool for the analysis of 
future measurements, we also calculate the relevant 
observational parameters as a function of the inter- 
stellar plasma density. 

2. Formulation of the Problem 

The interface is characterized by three surfaces: 
the solar wind termination shock (TS), the heliopause, 
and the interstellar bow shock. The interstellar atom 
flow in the heliospheric interface must be described 
kinctically because the mean free path of the neutral 
atoms is of the order of the size of the heliospheric 
interface. In fact, Baranov et al. [1998] have shown 
kinetic and hydrodynamic models of the interstellar 
atom flow may give significant differences. In order to 
obtain the kinetic distribution function, the Boltzman 
equation must be solved: 

_ dfa(r, !?h) , F dfn(r,wn) 

WH 1 WZi = 

or bih owh 

-/h {f, wh) J | «fe - w p | a™ f p (r,w p ) dw p (2) 

+/ p (f, wh) J \w n , - whW™ fp(r, wh* )dwE_' 

+ ^impact )M(r, life) 

Here /h(^, wh) is distribution function of the H atoms, 
/ p (r, w p ) is the local distribution function of the pro- 
tons which is assumed to be Maxwellian, w p and % 
are the individual proton and H atom velocities, re- 
spectively. Here cr^P is the charge exchange cross 
section of an H atom with a proton, 0i is the pho- 
toionization rate, tor is the mass of atom, /3i mpac t is 
the electron impact ionization rate, and F is the sum 
of solar gravitational force and solar radiation pres- 
sure force. 

Equation (2) takes into account the following pro- 
cesses. 

l.The resonance charge exchange process: 

H + H + <-> H + + H 

with charge exchange cross section [Maker and Tins- 
ley, 1977] a™ = 1CT 14 (1.64 - 0.0695 In V) 2 , cm 2 . 
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Here V is the relative velocity measured in centime- 
ters per second. 

2. The photoionization process: The photoioniza- 
tion rate is 

/?f=/3 p V(7) 2 ,^V = 8.8-10- 8 s- 1 

where 0^ h E is the photoionization rate at the Earth's 
orbit and r e is 1 AU. 

3. The electron impact ionization process: The 
rate is given by [Lotz, 1967] 

fc = £y^«VA W A) (3) 

-b-e'-^EkiX + c)} 

Here A = Ei(X) = f™ ^£&dt. The values 

Pi,a,b,c are equal to Pi = 13.6 eV, a = 4. • 10~ 14 
cm 2 cV 2 , b=0.6, c=0.56. 

4. The solar gravitation (F g ) and solar radiation 
pressure ( F r )processes: 

F = F g - F r = (1 - fi)F g , n = F r /F g 

We have used the value /U = 0.8. Most of our results 
are independent of the chosen value, since F has a 
nonnegligible influence only within a few AU from 
the Sun. 

The difficulty with modeling the H atom flow lies 
in the necessity to take into account the mutual influ- 
ence of the atomic and plasma components [Baranov 
and Malama, 1993, 1995, 1996] and to solve the ki- 
netic equation (2) together with hydrodynamic equa- 
tions for the plasma component. To calculate the 
H atom flow, we used the axisymmetric model and 
the method developed by Malama [1991] and Baranov 
and Malama [1993, 1995, 1996]. The boundary condi- 
tions for the proton density, the bulk velocity and the 
Mach number of the solar wind at the Earth's orbit 
are taken as n Pv E = 7.00 cm' 3 , V# = 450 km s _1 , 
M B =10. 

In the unperturbed LIC, we use Vlic = 25 km 
s _1 and Tlic = 5600 K for all sets of model param- 
eters. These values are close to the most recent de- 
terminations of interstellar He parameters obtained 
by Witte et al. [1996] with the GAS instrument on 
Ulysses. These authors give an interstellar helium ve- 
locity 24.6 ± 1.1 km s _1 and a helium temperature 
of 5800 ± 700 K. The H atom number density is kept 
fixed at nn,Lic = 0.2 cm" 3 , as discussed in the intro- 
duction. This value corresponds to the mean value 



given by Gloeckler et al. [1997]. For our calculations 
we have chosen the following values of proton number 
density: n PiL i C = 0.3,0.2,0.1,0.07,0.04 cm- 3 . 

3. Results of Modeling Calculations 

Using the heliospheric interface model with the so- 
lar wind and LIC parameters described above, we 
have calculated the structure of the heliospheric inter- 
face (positions and shapes of the TS, HP, and BS) and 
the distributions of plasma and neutral components. 
The distributions of the plasma and of the different H 
atom populations as well as the influence of the differ- 
ent physical processes have been discussed by Bara- 
nov and Malama [1993, 1995, 1996] and Baranov et al. 
[1998]. Here we present and discuss only selected re- 
sults of our calculations, which will be useful for anal- 
yses of interstellar proton number density in the sec- 
tion 4. Figure 2 shows the BS, HP,and TS distances 
to the Sun in the upwind direction as a function of 
the interstellar proton number density. It can be seen 
from the figure that the BS has the largest response to 
the interstellar proton density variations. As a matter 
of fact, the distance to the BS in the upwind direction 
decreases from 360 AU (for n Pi Lic = 0.04 cm -3 ) to 
180 AU (for n p , lic = 0.3 cm -3 ), whereas the distance 
to the heliopause varies from 185 to 110 AU. Thus the 
heliospheric interface region (the region between the 
BS and TS) becomes narrower while the interstellar 
proton number density is increased. The distance to 
the termination shock is NOT very sensitive to the 
proton density (between sa 100 and w 70 AU). The 
inferred range is compatible with the shock location 
deduced from the radial gradients of the ACR energy 
spectra [e.g. Cummings and Stone, 1996]. 

400 V J 




F 1 

0.0 0.1 0.2 0.3 

n P ,Lic cm ~ 3 

Figure 2. Positions of the bow shock, the heliopause and the 
termination shock in upwind direction as a function of interstellar 
proton number density. 



5 



0.50 




100 200 300 400 

Distance, AU 

Figure 3. The number density of interstellar atoms as a function 
of heliosphcric distance for different values of interstellar proton 
number density. 



cm 3 , the values are Th,ts = 14000 K and Vh.ts = 
17.0 km s _1 . 




0.0 0.1 0.2 0.3 



n p,LIC Cm 



Figure 3 shows the number density of interstellar 
H atoms as a function of the heliocentric distance in 
the upwind direction. As the proton number density 
increases, the hydrogen wall between the BS and the 
HP becomes "higher." The filtration factor (defined 
in equation (f )) decreases. Indeed, an increase of the 
proton number density in the LIC leads to increase 
of the proton number density between the BS and 
the HP and to an increase of the secondary H atoms 
number density resulting from the charge- exchange 
between primary interstellar H atoms and decelerated 
protons. This population of secondary H atoms has a 
smaller bulk velocity and a higher temperature than 
the primary interstellar H atom population. It is the 
increase of the number density of the secondary H 
atoms which reinforces the H wall and decreases the 
filtration factor. Figure 3 also shows that filtration 
occurs mainly in the region between the BS and the 
HP. 

Table 1 shows how the number density of primary 
interstellar and secondary H atoms at the termination 
shock in the upwind direction changes as a function of 
the density n Pi Lic- For n p ,Li C = 0.04 cm~ 3 , primary 
interstellar H atoms represent about 50% of the H 
atoms entering the heliospheric interface, while for 
?ip,Lic = 0.2 cm~ 3 , they represent only 6% of the 
total. Since the secondary H atoms have a smaller 
velocity and a higher temperature than the primary 
interstellar atoms [cf., Baranov et at, 1998], the bulk 
velocity and the temperature of the mixed H atom 
gas vary with the interstellar proton number density 
(Table 1). For example, for n Pi Lic = 0.04 cm~ 3 , the 
temperature and the velocity are Th,ts = 10500 K 
and Vh,ts = 22.5 km s -1 , while for n P; Lic = 0.3 
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0.0 0.1 0.2 0.3 




0.0 0.1 0.2 0.3 



n P ,uc c m 

Figure 4. (a) Number density, (b) the velocity, and (c) the tem- 
perature of intertcllar atoms at the termination shock in upwind di - 
rection as functions of interstellar proton number density. 

Figures 4 shows the interstellar atom density, ve- 
locity, and temperature at the termination shock in 
the upwind direction as a function of n Pj Lic- It can be 
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seen from the figure that the interstellar atom num- 
ber density and the filtration factor decrease rapidly 
when the proton number density is increased from 
to 0.04 cm -3 and much less rapidly for higher n p L ic 
values. As a consequence, small uncertainties in the 
atom density at the TS will give us small uncertain- 
ties in n Pi Lic for n p,Lic < 0.04 cm" 3 ("low density 
case") and relatively large ones for, say, n Pi Lic > 0.07 
cm" 3 ("high density case"). The same situation oc- 
curs for the temperature (Figure 4c), showing that 
its determination is an excellent diagnostic in the 
low density case. Contrary the filtration factor and 
temperature, the velocity dependence on n Pi Lic (Fig- 
ure 4b) remains about the same for the full range 
cm" 3 < n p lic < 0.3 cm" 3 , showing that knowledge 
of Vh,ts is very helpful for determining n Pj Lic even 
in the "high density case" (as long as n Pi Lic < 0.2 
cm -3 ). 

Figure 5 shows the plasma frequency in the inter- 
stellar plasma and at the bow shock on the upwind 
axis as a function of the proton density. These fre- 
quencies may be important if, as it has been suggested 
by Gurnett et al. [1993], the 1.8 kHz emission cutoff 
corresponds to the interstellar plasma frequency, and 
if, as suggested by Grzedzielski and Lallement [1996], 
the 2 kHz band was emitted in the compression region 
ahead of the bow shock. Also plotted is the frequency 
corresponding to the maximum plasma density along 
the upwind and crosswind axis (the "pile-up" region), 
which according to the present model occurs between 
the BS and the HP. This may be an important param- 
eter too, since, according to the above scenario this 
is the "obstacle" the 2 kHz signal has to overcome to 
be able to enter the heliosphere. 

Recently, Linsky and Wood [1996] have convinc- 
ingly shown that the excess of neutral hydrogen ab- 
sorption seen in the spectrum of the star a Centauri 
had indeed its origin in the heated and decelerated 
gas from the so-called H wall corresponding to the 
compressed region between the BS and the HP. While 
being an important discovery, simulations show it is 
hard to distinguish from these observations between 
the H walls built up for n Pi Lic = 0.04 or 0.2 cm" 3 . As 
a matter of fact , if n Pi Lic increases, the gas is more 
heated and compressed, but the thickness of the H 
wall is reduced. This is why we have not included the 
H wall absorption in this parametric study. A some- 
what larger difference does exist between the "super- 
sonic" and "subsonic" cases, as calculated by Gay- 
ley et al. Here we consider the supersonic case only. 




Figure 5. The plasma fro quencics in the interstellar plasma 
(solid curve 1), at the bow shock (curves 2 and 3), and in the ma- 
ximum plasma density region between BS and HP (curves 4 and 5) 
as functions of interstellar proton number density. Dotted curves 
3 and 5 correspond to upwind, dashed curves 2 and 4 correspond 
to crosswind. 

4. The Influence of the Interstellar 
Proton Number Density 

4.1. The Neutral Hydrogen Density in the 
Inner Heliosphere 

The pick-up ions measurements provide a determi- 
nation of the neutral H flux which has the advantage 
of being dependent on the solar ionization processes 
only. The situation is even better since the solar wind 
measured at the same time as each set of pick-ups has 
indeed been the main ionizing agent of these partic- 
ular pick-ups. At variance with the pick-ups, the use 
of the backscattered Ly a glow as a neutral H density 
diagnostic suffers from uncertainties on photon radia- 
tive transfer effects and on radiation pressure and ion- 
ization rate measurements [see, e.g., Quemerais et at, 
1994]. 

As mentioned above, Gloeckler et al. [1997] have 
obtained from the SWICS Ulysses pick-up ion mea- 
surements for the neutral hydrogen number density 
at the termination shock 

0.09 cm" 3 < rijj.TS < 0.14 cm" 3 (4) 

using the classical "hot" model. Because the "hot" 
model does not take into account any heliosphcric fil- 
tering at the interface, it is reasonable to assume that 
this measured value corresponds to interstellar atom 
number density at the TS. Comparing this with the 
results of our numerical calculations (Figure 4a), we 
obtain 

0.02 cm" 3 < n Pi Lic < 0.1 cm" 3 (5) 
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Table 1. Interstellar Atoms Number Density at the Termination Shock (tih,lic = 0.2 cm~ 3 ) 



«p,lic, cm 3 n prim ,TS, cm 3 n SCCjT s, cm 3 n H ,TS, cm 3 Vb,ts, km s 1 T h ,ts, K 



0.3 0.07 17 14000 

0.2 0.0045 0.075 0.0795 18 13500 

0.1 0.02 0.07 0.09 20 12500 

0.07 0.03 0.065 0.095 21 12000 

0.04 0.055 0.05 0.105 22.5 10500 



Table 2. Intervals of Possible Interstellar Proton Number Densities 



Type of Hcliosphcric Interface Diagnostics 



Range of Interstellar Proton Number Density 



SWICS/Ulysses pick-up ion 
0.09 cm" 3 < n H ,TS < 0.14 cm" 3 
Gloeckler et at, [1997] 

Ly- a, intensity 
0.11 cur 3 < u h ,ts < 0.17 cm" 3 
Quemerais et al., [1994] 

Ly-a, Doppler shift 

18 km s" 1 < Vh^ts < 21 km s" 1 

Bertaux et al. [1985], Lallement et al, 

[1996], Clarke et al. [1998] 
Voyager kHz emission (events) 

110 AU < Rav < 160 AU 

Gurnett and Kurth [1996] 

Voyager kHz emission (cutoff) 
1.8 kHz 

Gurnett et al. [1993], Grzedzielski and 
Lallement [1996] 



0.02 cm 3 < n P! Li C < 0.1 cm" 



"p.lic < 0.04 cm 3 or 

«p,lic < 0.07 cm" 3 (for n H ,Lic = 0.23 cm" 3 ) 
0.07 cm" 3 < n P! Li C < 0.2 cm" 3 



0.08 cm" 3 < n p ,Lic < 0.22 cm" 



"p.lic = 0.04 cm" 
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See, also, table 2. The lower value of the neutral den- 
sity interval falls within the insensitive part ("high 
density case") of the function nH,Ts( n p,Lic) in fig- 
ure 4a, and this explains the large range resulting 
for ro Pi Lio This is something only estimates as those 
from Gloeckler et al. [1997] do not show and the 
present work makes visible. The mean Gloeckler et 
al. value njj,TS = 0.115 cm -3 (nH,Ts/^H,Lic = 0.575) 
corresponds here to n p .Lic = 0.03 cm" 3 . This value 
is slightly lower than the one obtained by Gloeckler 
et al. In order to narrow the range for the interstellar 
proton density on the basis of pick-up measurements, 
we need a larger precision on the neutral H density. 
As we mentioned above, the backscattered solar Ly a 
intensity can also provide an estimate of the interstel- 
lar atom number density. A compilation and reinter- 
prctation of many previous determinations of the H 
atom density by Ly a measurements has been done 
by Quemerais et al. [1994]. The resulting density is 
in the range of 0.11 cm" 3 - 0.17 cm" 3 , corresponding 
to 

n p , L ic < 0.04 cm" 3 (6) 

The range obtained by Quemerais et al. includes the 
mean value of the n p ,Lic determination by Gloeckler 
et al. 

4.2. Velocity and Temperature 

Figure 4 shows that the temperature and the veloc- 
ity of the H atoms in the heliosphere depend rather 
strongly on the LIC proton number density. These 
model predictions can be compared with estimations 
obtained from measurements of the backscattered so- 
lar Ly a radiation. Velocities are deduced from 
Dopplcr shifts, while temperatures are deduced from 
linewidths with the help of flow modeling. However, 
one has to overcome two difficulties: 
(1) The velocity determination suffers from uncertain- 
ties in the radiation pressure. Depending on the bal- 
ance between radiation pressure and gravitation, the 
velocity of the gas close to the Sun increases (p< < 1 
valid for low activity) or decreases (pi > 1 valid for 
high activity). Therefore an uncertainty in the bal- 
ance parameter (i introduces an error in the H velocity 
far from the Sun. (2) The temperature determination 
suffers from uncertainties on the radiative transfer ef- 
fects, which broaden the lines in a manner that is not 
yet satisfyingly represented. 

Measurements of the interplanetary Ly a emission 
line profile obtained with a hydrogen absorption cell 
on board Prognoz 5/6 [Lallement et at, 1984; Bertaux 
et al., 1985] have yielded a good estimate of the bulk 



velocity (20 ± 1 km s" 1 ). In this case, fi was de- 
termined by adjustment of the model. However, the 
"hot" model never reproduced the data in all direc- 
tions simultaneously. More recently, Ly a spectral 
measurements [Lallement et al., 1996, Clarke et at, 
1998] confirmed that the inflow speed of H atoms far 
from the Sun (at the TS) is within the range of 18 - 
21 km s" 1 using estimates of /i based on solar Ly a 
measurements corresponding to the periods of obser- 
vation. 

According to the results of our calculations (Figure 
4b), the velocity interval 18 km s" 1 < Vh,ts < 21 km 
s" 1 corresponds to the range 

0.07 cm" 3 < n PiL ic < 0.2 cm" 3 (7) 

The analysis of the Prognoz 5/6 H cell data also 
provided the "line profile" temperature (8000 ± 1000 
K) of the H atoms in the inner heliosphere. This 
is about 1500-2500 K above the helium temperature. 
However, it is important to note that only line-profiles 
observed in the direction perpendicular to the main 
flow could be used and traced back through modeling 
to the temperature "at infinity," i.e., before the inter- 
action with the Sun (but still inside the heliospheric 
interface). In addition to the uncertainties on the ef- 
fect of radiative transfer on line profiles, this determi- 
nation may have been slightly biased by the assump- 
tion of a given temperature in the absorption cell and 
the temperature slightly underestimated. Now, the 
Hubble Space Telescope Ly a spectral observations 
mentioned above have shown that the line profiles on 
the upwind and downwind side are much larger than 
what predict classical models for a gas initially at the 
same temperature as helium. Inferred line-of-sight 
temperatures are as high as 15,000-20,000 K [Clarke 
et al, 1998]. This is in favor of a nonncgligiblc heating 
of hydrogen at the interface. Still, the role of radia- 
tive transfer has to be assessed before one can derive 
a confidence interval for the kinetic temperature. 

4.3. Heliospheric Radio Emissions 

The radio emission detected by Voyager 1 and 2 
spacecraft is another type of heliospheric interface di- 
agnostics. It is believed that the emission region is 
connected with the vicinity of the heliopause. Major 
events of kilohertz emission observed in 1983-1984 and 
in 1992-1993 were associated with intense solar wind 
solar events in 1982 and 1991. The propagation delay 
in both cases is about 400 days. Using measurements 
of the propagation speed of the interplanetary shocks 
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produced by these events and the time delay of the on- 
set of the radio burst, Gurnett and Kurth [1996] could 
estimate the distance to the heliopause and found that 
it ranges from 110 to 160 AU. The comparison with 
our calculations (Figure 2) shows that this range for 
the HP location (for the upwind side) corresponds to 
an interstellar proton number density 

0.08 cm" 3 < rc Pi Lic < 0.22 cm" 3 (8) 

The second feature is the lower frequency (1.8-2.1 
kHz) emission band and in particular the well-defined 
"cutoff" of this emission at 1.8 kHz. This sharp "cut- 
off" could be related to the LIC density, which is the 
unique constant parameter for the whole interface. 
The plasma frequency at 1.8 kHz implies an interstel- 
lar electron density of 0.04 cm" 3 . As explained in the 
previous section, any radiation emitted in the inter- 
stellar medium is, in principle, prevented to enter the 
interface due to a maximum compression region char- 
acterized by a density we have represented in Figure 
5. Despite these difficulties, it remains that there is 
no other explanation for the cutoff, and the particu- 
lar value n p ,Lic = 0.04 cm" 3 has a high probability 
of being the true circumsolar interstellar density. 

4.4. A Synthesis? 

It can be seen from equations (6) and (7) that 
there is already a small discrepancy between the in- 
terstellar plasma densities obtained from Ly a inten- 
sity measurements and the Ly a profile (Doppler's 
shift) measurements. At the same time, both inter- 
vals have an intersection with the interval in equation 
(5) derived from pick-up ions measurements. This 
discrepancy disappears if both the interstellar pro- 
ton and neutral number densities are higher than the 
values we have assumed. In this case, it is possible 
to have simultaneously a smaller bulk velocity and 
a higher neutral H density in the hcliosphere. How- 
ever, there is a model-independent limit on the LIC 
H atom density from the relative H/He abundance 
«h,lic < 0.23 cm" 3 (see section 1). This implies that 
only a small increase of tih.lic (and then of n P) Lic) 
is relevant. For the maximum value nn,Lic = 0.23 
cm" 3 and n p ,Lic = 0.07 cm" 3 (in agreement with 
equation (5)), we find nn.TS = 0.11 cm" 3 , assuming 
the same filtration factor as for Uh.lic = 0.20 cm" 3 
and n p xic = 0.07 cm" 3 . This value corresponds to 
lower limit derived from Ly a intensity measurements. 

Up to now, we can conclude that the pair of values 
«h,lic = 0.23 cm" 3 and n Pi Lic = 0.07 cm" 3 is in 
agreement with both observations of pick-up ions and 



Ly a radiation. Now, we note that n Pi Lic = 0.07 
cm" 3 is very close to the lower limit in equation (8). 
Considering that the source of the 3 kHz emission may 
not be exactly at the HP but in the region between the 
TS and the HP could help to reconcile the two results. 
This could be due to, for example, the influence of the 
electron impact ionization on the plasma flow in the 
region between the TS and the HP. As a matter of 
fact, this process leads to a strong plasma density 
gradient between the HP and the TS. 

However, in any case, the above values are not com- 
patible with the 1.8 kHz cutoff, if it corresponds to 
the interstellar emission. If we assume now that this 
is really the case, then we have to explain why equa- 
tions (6) and (8) are not justified. 

We can give two possible explanations for the dis- 
crepancy with equation (8). First, it is possible that 
the source of the emission is not exactly at the HP 
but in the region between the TS and the HP, as it 
was mentioned above. A second possible explanation 
is that there is some additional pressure in the inter- 
stellar medium. It may be a magnetic field pressure 
[Myasnikov, 1997; Linde et al., 1998; Pogorelov and 
Matsuda, 1998] or a low energy cosmic ray pressure 
[Izmodenov, 1997; Myasnikov et al., 1997]. Owing to 
this additional pressure, the HP would be closer to 
the Sun as compared with what the present model 
implies for n Pi Lic as low as 0.04 cm -3 . 

The calculation of the interstellar magnetic field 
(IMF) strength required to push the HP as close as 
110-160 AU requires a full MHD model coupled to a 
neutral flow model. Such a model is not yet available. 
The magnetic field is taken into account in some gas- 
dynamical models by modifying the equation of state 
for the plasma [e.g., Gayley et al, 1997], which is 
appropriate in some specific conditions of orientation 
and Mach numbers. MHD models have been built but 
without inclusion of the coupling to the neutral flow. 
Hereafter, we estimate the IMF strength needed to 
reconcile Hh.lic = 0.04 cm -3 , connected with the 1.8 
kHz cutoff, and HP distance measurements by Gur- 
nett and Kurth [1996]. We use the following formula 
deduced from our calculations and Figure 2 in the 
absence of magnetic and cosmic rays pressure: 

i? H p -24(3.6 -P L ic -10 12 ) + 110 (9) 

Here Php is the distance to the heliopause (upwind di- 
rection) in astronomical units, Plic is the interstellar 
pressure in dyn cm" 2 deduced from our parameters. 
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We then replace Plic by the more general term 

-Plic = « p ,lic • (2fc b T L i C + toh^lic) ( 10 ) 
B 2 

+ a ^T^ + P GCR,LIC 
87T 

where n p .LiCj ^Licand Vlic are interstellar proton 
number density, temperature, and velocity, respec- 
tively, -Blic is the interstellar magnetic field strength, 
a is an amplification factor (a « — 4) determined by 
the angle between the interstellar magnetic field and 
velocity vector [Holtzer, 1989; see, also, Frisch, 1993], 
and Pgcr,lic is galactic cosmic ray pressure. 

By doing so, we assume that an additional non- 
gasdynamical pressure acts on the shape of the he- 
liopause on the same way as a gasdynamical (proton 
and electron) pressure. The influence of interstellar 
atoms is included in the numerical coefficients of equa- 
tion (9), which implies this formula is valid only if 
«h,lic is of the order of 0.20 cm~ 3 . 

It is interesting to note that if the direction of the 
interstellar magnetic field is the same as the wind di- 
rection (in principle the unique possibility for a two- 
dimensioned modeling), the heliocentric distance of 
the heliopause increases as a result of the magnetic 
field tension [Baranov and Zaitsev, 1995] along the 
wind axis, whereas this distance decreases in the wing 
due to magnetic field pressure. In this case, the for- 
mula (9) does not work. However, the parallel IMF is 
probably unlikely because, to reconcile data, we need 
the heliopause closer to the Sun. 

Our formula (9) probably gives better estimates of 
the HP distance variation as a function of the inter- 
stellar pressure than estimates based on balance be- 
tween the solar and interstellar pressures [see, e.g., 
Holtzer, 1989], because we take into account the 
plasma compressibility and the influence of interstel- 
lar neutrals on the plasma structure. Actually, owing 
to these effects, the heliospheric interface is a kind 
of damping region for any changes in the interstellar 
pressure. Using balance pressure only, one overesti- 
mates the variation of the HP distance due to inter- 
stellar pressure. 

If we assume that n p ,Lic = 0.04 cm~ 3 and that 
the GCR pressure is Pgcr = 0.3 • 10~ 12 dyn cm~ 2 
[Holtzer, 1989], we find from equation (9) that the in- 
terval 110 AU < i? H p < 160 AU [Gurnett and Kurth, 
1996] corresponds to 

B 2 

0.737- 10- 12 dyncm- 2 < a— < 2.82-10~ 12 dyn cm" 2 

87T 



For a— A (perpendicular magnetic field) [Holtzer, 1989], 

2.1 /iG < B < 4/zG. 

This interval is in agreement with the current esti- 
mates of the interstellar magnetic field strength [e.g., 
Frisch, 1995]. However, it is necessary to note that 
if indeed the IMF strength is above 2.1 /zG, then a 
full "subalfvenic" model is required. However, in the 
subalfvcnic case, there is no bow shock and then, to 
date, no explanation for the 2 kHz radio emission. 

It remains that despite the addition of such a pres- 
sure, the discrepancy with the bulk velocity measure- 
ments equation (6) still remains. If n p xic = 0.04 
cm~ 3 , then the bulk velocity is Verts = 22.5 km s _1 
(23 km s _1 for an interstellar velocity of 25.5 km s _1 ). 
As we already discussed, the Doppler shift measure- 
ments are sensitive to the Ly a radiation pressure. 
The above value of the bulk velocity implies that the 
radiation pressure is above what has been inferred 
from the H cell data (there is less deceleration in- 
duced by the interface and more due to the radiation 
pressure). 

All the above estimates have been done by im- 
plicitely assuming that the additional pressure due 
to the fraction of neutral gas coupled to the plasma 
remains of the same order with and without interstel- 
lar magnetic field. Of course, both the intensity and 
the direction of the magnetic field change the plasma 
pressure between the bow shock and the heliopause, 
but they also change the thickness of this region. An 
increase of the plasma compression in all cases cor- 
responds to a decrease of the thickness; That is, the 
two phenomena tend to compensate, as it can be seen 
in the results of Linde et al. [1998] and Baranov and 
Zaitsev [1995], for direction and intensity changes, re- 
spectively. The filtering depends on the product of 
plasma density and thickness, and, as a consequence, 
it should not change dramatically. More accurate 
computations on the basis of MHD models are needed 
to quantify such changes. 

5. Conclusions 

We have performed a parametric study which shows 
how sensitive and compatible are the various types of 
diagnostics of the interstellar plasma density ™ Pj liCj 
i.e., the interstellar neutrals temperature, number 
density and velocity, and the radio emissions time de- 
lays and frequency ranges. For the neutrals, there are 
two regimes: For low values (n Pj Lic < 0.05 cm -3 ), 
the most sensitive parameters are the neutral H den- 
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sity and temperature in the inner hcliosphcrc, while 
for higher values, the H bulk velocity only remains 
sensitive. 

In the light of this study, we have discussed the 
observational results. Our main conclusion is that it 
is impossible to reconcile the results obtained from all 
types of data as they stand now. There is a need for 
some modifications of the interpretations or the con- 
fidence intervals. Two types of solutions (which are 
mutually exclusive) seem to be favored: (1) It is pos- 
sible to reconcile the pick-up ions and Ly a measure- 
ments with the radio emission time delays if a small 
additional interstellar (magnetic or low energy cosmic 
ray) pressure is added to the main plasma pressure. 
In this case, n Pi Lic = 0.07 cm~ 3 and Uh.lic = 0.23 
cm~ 3 is the favored pair of interstellar densities. How- 
ever, in this case, the low frequency cutoff at 1.8 
kHz of course cannot be connected to the interstel- 
lar plasma density, and one has to search for another 
explanation. (2) The low frequency cutoff at 1.8 kHz 
is connected to the interstellar plasma density, i.e., 
^p,Lic = 0.04 cm~ 3 . In this case, the bulk velocity 
deduced from Ly a spectral measurement is under- 
estimated by about 30-50% (the deceleration is by 3 
km s _1 instead of 5-6 km s _1 ). Model limitations 
(as the use of a stationary classical hot model to de- 
rive the bulk velocity [Rucinski and Bzowski, 1996]) 
or the influence of a strong solar Ly a radiation pres- 
sure may play a role. However, in this case, there is a 
need for a significant additional interstellar pressure 
as compared with case (1). If the source of this extra- 
pressure term is a perpendicular magnetic field, its 
strength should be in the interval 2.1 /xG < B < 4 /j,G, 
a value in agreement with local IMF estimates. How 
such an additional field will modify our conclusions on 
the interstellar plasma density is still an open ques- 
tion. However, we do not expect substantial changes, 
as we have discussed in the section 4. 

A need for an additional pressure is in agreement 
with the conclusions of Gayley et al. [1997] from the 
analysis of the H wall absorption toward alpha Cen- 
tauri. However, it remains that since the best model 
of these authors corresponds to a neutral H density of 
0.025 cm~ 3 in the inner heliosphere, at least 4 times 
smaller than the density derived from the pick-up 
ions, additional calculations for more realistic den- 
sities are still needed. 

New observations of the heliospheric gas and ions 
are expected within the next years. In particular, 
the SWAN instrument on board SOHO has gathered 
a considerable amount of Ly a data. Their analy- 



sis is in progress and will provide extremely precise 
measurements of line-of-sight temperature and bulk 
velocities of atomic H in all directions, as their varia- 
tions change year after year. This should help to dis- 
entangle the differences between the interpretations 
discussed here. In parallel, there is a crucial need for 
MHD model developments and for an unambiguous 
interpretation of the 1.8 kHz cutoff. 
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